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We study an ultracold atom-cavity coupling system, which had been implemented in experiment 
to display weak light nonlinearity [S. Gupta et at, Phys. Rev. Lett. 99, 213601 (2007)]. The 
model is described by a non-interacting Bose-Einstein condensate contained in a Fabry-Perot op- 
tical resonator, in which two incommensurate standing-wave modes are excited and thus form a 
quasiperiodic optical lattice potential for the atoms. Special emphasis are paid to the variation of 
atomic wavefunction induced by the cavity light field. We show that bistability between the atomic 
localized and extended states can be generated under appropriate conditions. 
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I. INTRODUCTION 

We consider the following model depicted in Fig. [T] 
A scalar Bose-Einstein condensate (BEC) with atomic 
number N is confined in a high-finesse Fabry-Perot cav- 
ity along the cavity axis in the x-direction, in which two 
standing-wave modes are excited. The atomic motion 
in the transverse direction is restricted to the ground- 
state of a strong harmonic trapping potential. At zero 
temperature the condensate can be described by the 
wavefunction ip(x,t). The interaction between the two 
cavity modes and atoms trapped inside is of dispersive 
nature. One cavity mode is relatively strong and ac- 
tively locked through cavity feedback, thus providing a 
time-independent primary optical lattice potential for the 
atoms. We call this mode as the trap mode, character- 
ized by the wavenumber kr, and which is not affected by 
the atomic dynamics. The second cavity mode with fre- 
quency uj c and wavevector k p = fikr, is relatively weak 
and is influenced by the atomic dynamics. Thus this 
mode can serve as a probe, which is driven by a coher- 
ent laser field with frequency ui p and amplitude rj. We 
will therefore call this mode as the probe mode. When 
/3 is an irrational number, the two cavity modes are in- 
commensurate with each other, and together they form 
a quasiperiodic lattice for the atoms. 




FIG. 1: (Color online) Schematic diagram showing the system 
under consideration. 

Such a model has already been realized in experiment, 



with the emphasis on the cavity electrodynamics [l[ and 
cavity optomechanics Q (which assumes that the atomic 
deviation from its equilibrium position is small). Due to 
the incommensurance of the trap mode and probe mode, 
atoms in each well of the primary lattice is subject to a 
different optical force induced by the probe mode. The 
collective effect in turn induces a Kerr-type nonlinear- 
ity (proportional to the intensity of the probe mode) to 
the probe-pump detuning. The strong nonlinearity thus 
allows for optical bistability at photon numbers below 
unity. 

On the other hand, it is well-known that the wave- 
function describing particle motion in a periodic poten- 
tial with disorder can generally be grouped into extended 
states and localized states. The transition between them 
was originally studied by Anderson in the system of non- 
interacting electrons in a crystal with impurities Q . An- 
derson localization describes the localization of waves in 
disordered media, which is considered as not only a fun- 
damental phenomenon in condensed-matter physics but 
also a universal phenomenon in many areas of physics, 
such as light propagation in photonic crystal [4j. The 
combination of ultracold atoms and optical potentials 
have provided a powerful playground for the study of 
disorder-related phenomena [5(- Specifically, in the case 
of an optical lattice a pseudorandom potential can be 
generated by superimposing two standing waves of in- 
commensurate wavelengths. The localization of atomic 
matter wave in such a quasiperiodic optical lattice have 
been experimen tally observed @, 0] and extensively stud- 
ied in theory |8l-[L5j. 

In the localized regime, atoms accumulate in one or 
several lattice sites, thus strongly deviate from its equi- 
librium position. Some questions then naturally arise: 
Whether the localization of the atomic wavefunction can 
be observed in the atom-cavity coupling system described 
above? What is the effect of cavity on localization? Will 
optical bistability still take place in the localized regime? 

Although the study of interplay between the atomic ex- 
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ternal or internal degrees of freedom and the cavity light 
field have already attracted a lot of theoretical jl7H28| 
and experimental interest [29l - l3lj | , the new interest in this 
model lies in how the localization transition of the atoms 
gets modified inside the optical resonator. The quasi- 
disorder induced by the probe mode now depends on the 
atomic distribution, while the properties of the atomic 
wavefunction is intimately related to disorder. The prob- 
lem is then highly nonlinear. Through this work, we ex- 
pect to shed some light on the problems related to the 
interplay between nonlinearity and disorder. 

The rest of the paper is organized as follows. Sec II 
introduces the theoretical model, we show that the sys- 
tem can be mapped into the Aubry- Andre model in the 
tight-binding limit. Sec III is devoted to the discussion 
of a special case, in which the physical properties of the 
system can be captured by a cavity-dressed asymmetric 
double-well system. Both the equilibrium properties and 
dynamics are discussed. The more general case is con- 
sidered in Sec IV, in which we present numerical results 
and display the bistability between the atomic extended 
and localized states. Finally we conclude in Sec V. 



II. MODEL 

In order to concentrate on the nonlinearity induced by 
the cavity light field, we assume that the atomic s-wave 
scattering leng th is tuned to zero by means of Feshbach 
resonance [32j. With the energy measured in units of 
tvjjR = h 2 kj j /2m (wji ~ 2tt x 5 kHz for 87 Rb atoms) and 
length scaled in units of fc^ 1 , the Gross-Pitaevskii (GP) 
equation describing the evolution of the condensate wave- 
function and the corresponding mean-field equations of 
motion for the cavity light field read 



dx 2 



„ + U t sin 2 (x) 



+ \a(t)\ 2 U sin 2 {fix + (/>) ip(x,t), (la) 

—a (t) = \iS c - iNU (sin 2 (/3x + 4>)) - n] a ft) + r\. 
at 

(lb) 

Here S c = uj c — uj p is the cavity-pump detuning, U t (n) 
characterize the strength of atomic coupling with the trap 
mode and probe mode, respectively. <fi is the phase dif- 
ference between the two cavity modes, a represent the 
amplitude for the probe mode with the decay rate k. 
For the simplicity of the following discussion, we assume 
U t (o) > an d <j> = 0. The bracket in Eq. (flbf is defined 
as 

(/(x)) = / dx|V(x)| 2 /(x). 



In the tight-binding limit, the atomic wavefunction 
tp (x) can be expanded over the lowest band Wannier ba- 
sis Wj (x) = w (x — j7r) defined by the primary lattice, 



i.e., ip (x) = J2j c j w j ( x ) w hh ^ - \cj\ = 1. Retaining 
only the coupling between neighbouring Wannier states 
and the onsite contribution of the probe mode, after ne- 
glecting the constant terms, Eq. (fTa|) can be rewritten 
as 



icj = — J (cj+i + Cj-_i) — A cos (2nj3j) Cj, (2) 



with 



J = — / dx Wj + i (x) 



-— 2+ U t ^{x) 



Wj (x) 



A= H_Uo j dxcQs p0 x )\ w ( x ) | 2 = x U \a\ 2 . 



In estimating the value of J and A, we take the Gaussian 
approximation for the wannier functions, which is valid 
for a sufficiently deep primary lattice, and from which we 
have 



J = ^Uf fi exp(-2y/U t ) • 

\/7T 



A 



H 2 u Q 



exp(-/3 2 / v / C/ t ) 



Taking the transformation 



(3a) 
(3b) 

(4) 



Eq. ([2j) can be cast into the form 



A 



(frn+l + fm-l) ~ 2 J COS (2tt/3j + k) f m . (5) 



Equation ([5]) has precisely the same form as Eq. ([2]) at 
A/ J = 2, which indicates the self-duality of the Aubry- 
Andre ( AA) model [33[ ■ Since Eq. (U]) represents the typi- 
cal discrete Fourier transform which transforms localized 
states into extended states and vice versa, A/ J = 2 is 
thus regarded as the transition point between the local- 
ized states and extended states. 

The linear case without atomic feedback on the cavity 
mode is exactly the same as those had been illustrated 
in Q , one can observe the AA localization by increasing 
the value of A/ J. The AA localization resembles the 
Anderson localization of random systems and have been 
experimentally observed with a non-interacting BEC in 
a bichromatic optical lattice @. In the following, we 
will investigate the nonlinear effect on AA localization 
brought out by the atom-cavity interaction. 



III. SUPERLATTICE 



= Cj . In 



Equation ([2J indicates the periodicity c J+ i 
order to gain some physical insights, we first consider 
the special case that 1//3 is an integer. Specifically, let 
us focus on the simple case where 1/(3 = 2. The cavity 
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FIG. 2: The superlattice potential formed by the cavity light field 
for /3 = 0.5, Ut = 10 and Uo \ce\ 2 = 1. The stripes indicate 
the first and second Bloch bands of the primary lattice formed by 
the trap mode. The pertubation introduced by the probe mode 
is indicated by the dashed line. The label 0(1) indicates the deep 
(shallow) site in a supercell. 



modes then form a superlattice potential (more specif- 
ically, a lattice of double-well potentials) for the con- 
densed atoms, as shown in Fig. [p Such a potential has 
been realized in the experiment [Iff. 

We have checked that, for the parameters considered 
here, the tight-binding assumption is valid and all the 
population are in the lowest Bloch band. Owning to 
the periodicity, the physics can be captured by one su- 
percell in the region of x £ [— A/4,3A/4], which repre- 
sents a cavity-mediated double- well system for the ultra- 
cold atoms. Similar systems have also been discussed in 
Ref. [26| . Due to the two independent cavity mode used 
here, the double-well potential is asymmetric, character- 
ized by a deep well and a shallow well 1. 



A. equilibrium properties: bistability between 
extended and localized states 

Following Eq. ((2J, the effective equations of motion 
inside one supercell can be written as 

icQ = — 2 Jc\ — Aco, 

ic\ = —2Jcq + Aci, (6) 

where we have used the periodic boundary condition 
c J+ 2 = Cj. This equation illustrates a non- interacting 
BEC inside a double-well potential with the tunneling 
coefficient 2 J and the trap asymmetry 2 A. Consider the 
following transformation according to Eq. Q 

CO = fo + fx, c\ = fo - fl, 



the equations of motion for fo and f\ reads 

ifo = -AA - 2J/ , 

i/i = -A/ + 2 J/i. (7) 

Eqs. (J7J) has precisely the same form as Eqs. © 
with the role of J and A interchanged, the two equa- 
tions are identical at A/ J = 2, which signal the crit- 
ical point. The physical implication of the ciritical 
point can be captured via the dimensionless parameters 

z — 2 d x I 7 /" ( x )\ 2 ~ 1 ~ l c o| 2 — I c i| 2 j which indicate 
the popultion imbalance between the two wells. Under 
the initial condition of cq (0) = 1, c\ (0) = (i.e., with all 
atoms initially prepared in the deep well), in the linear 
case where J and A are fixed external parameters, we 
have 

z(t) = l * sin 2 v / A 2 + 4J 2 i. (8) 

1 + (A/2 J) 

With the increase of A/2 J from to oo, Eq. (jHJ shows 
that the minimum value of z increases from —1 to 1, 
indicating that the atomic wavefunction become localized 
and are less likely to diffuse. At the critical value of 
A/ J = 2, the value of z oscillates between and 1, which 
means that at most half of the atoms can transport to 
the shallow well. In this sense we can group the atomic 
wavefunction into the localized state and extended state. 

In the nonlinear case, we will have to take the cavity 
feedback into account. According to Eq. (fTbf 

^-a (t) = [iA c + iNUqxz - k] a (t) + r), (9) 
at 

with A c — 6 C — NUq/2. By considering that the cavity 
decay rate k is typically much larger than the atomic 
oscillation frequency J, we adiabatically eliminate the 
probe field a from Eq. ([9]), replacing a in Eq. (|3b[) with 

OL= r- — -. (10) 

k - i [A c + NUaxz] 

Equations © can be rewritten in terms of z and the 
phase difference tp between the amplitudes Cq and C\ as 

z = 4J-\/l — z 2 sirup, 

<P = -4J^=L=cos^ + 2 x £/oM 2 . (11) 

The steady-state solution of the system can then be 
solved according to the above equations by takind the 
time derivative to be zero. An example is shown by 
the red dashed line in Fig. |3l which displays how the 
mean photon number of the probe mode changes with 
the increase of pumping amplitude tj/k. We can find 
three steady-state solutions in certain parameter region, 
two of them are dynamically stable while the third one 
is dynamically unstable, the unstable solution links the 
two stable ones, representing a typical dispersive optical 
bistability. 



FIG. 3: (Color online) The steady-state intensity of the probe 
mode |a ss | 2 versus r\j K using the two- mode treatment (red-dashed 
line) and the self-consistent imagninary-time propagation (black 
dotted line). The inset show (x)\ 2 (solid line) versus the optical 
potential (dashed line) inside one period for two different pump 
amplitude located at the two sides of the bistable region. The 
vertical black dashed line indicate the value of \a S s [ 2 corresponding 



to A/ J = 2. Parameters: 
U = 0.2uj r . Ut = Wujr. 
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The steady-state behavior can be understood as the 
interplay between the double-well asymmetry and the 
cavity-induced nonlinearity, based on a positive feedback 
mechanism: With the increase of the pumping amplitude 
rj for the probe mode, its intensity increases as well. It 
results in a larger A and thus enhance the double-well 
asymmetry felt by the atoms. As a result, the atoms 
tend to stay in the deep well. This will give rise to a 
larger population imbalance z and shift the cavity reso- 
nance, and will again lead to the increase of the probe 
mode intensity \a\ , as can be seen from Eq. (|10[) by 
setting A c < ~NU X- 

In order to understand the variation of atomic wave- 
function associated with the optical bistability, we obtain 
the system steady-state in a self-consistent manner by 
propagating Eq. (Ila|) in imaginary-time with the steady- 
state value of the probe field amplitude from Eq. (flbf . 
The numerical simulation is performed on one period 
A with periodic boundary condition. The result is dis- 
played as the black dotted line in Fig. [3j which shows 
good agreement with the two- mode result. The typical 
wavefunctions are shown in the insets. In the linear mode 
without cavity feedback, as the number of the probe pho- 
ton increases, the atomic system will experience a smooth 
crossover from the state that the atoms evenly populated 
between the two wells to that only populated in the deep 
well. Here in the cavity system, as the probe pumping 
rate rj/n is increased, the probe photon number will show 
discontinuous jumps. Accompanying with this jump in 
photon number, the atomic wavefunction also exhibits 
jumps between an extended state and a localized (self- 



FIG. 4: (Color online) (a) The steady-state intensity of the probe 
mode |a ss | 2 and (b) Inversion z versus cavity- pump detuning S c - 
The red-dashed is the result from the two-mode treatment while 
the black-dotted line results from the self-consistent imaginary time 
propagation. The parameters are the same as Fig. [3] except that 
77 = 120ujr. 

We should point out here that the dispersive optical 
bistability is not necessarily accompanied by the transi- 
tion of the atomic wavefunction across the critical point. 
An example is shown in Fig. 2J Here we fix the pump- 
ing amplitude 77 while varying the cavity-pump detuning 
5 C , as those had been done in experiment It is clear 
that the region for optical bistability is located below the 
critical point, meaning that the bistability is taken place 
between extended states. We anticipate that similar sit- 
uation appears in the work of Ref. 0, 0, in which the 
system operates with a weak probe light field, the atomic 
deviation from its equilibrium position is small and thus 
validates the mapping into cavity optomechanics. 



B. nonequilibrium dynamics 

In order to understand how the system dynamics are 
modified by the nonlinear interaction, we use two differ- 
ent numerical methods: (i) integrate the GP equation 
(jlap in a self-consistent manner by assuming that the 
probe field amplitude adiabatically follow the the vari- 
ation of the condensate wavefunction; (ii) use the two- 
mode model, and solve Eqs. (jlip with Eq. (TTU| . To 
make a comparison with the dynamics predicted in the 
linear case as represented by Eq. ([S]). we assume that 



5 



the atoms are initially localized in the deep well. The 
numerical results are shown in Fig. [5] 
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FIG. 5: (Color online) Evolution of population imbalance 2 and 
intensity of the probe mode \a\ 2 with time t for tj/k, — 0.5 (upper 
panel) and 1.5 (lower panel). The other parameters are the same 
as those in Fig. [3] The black solid line are the results calculated 
with numerical mithod (i) and red dashed is with method (ii). The 
dash dotted line is given by Eq. J5}. 



FIG. 6: (Color online) The pump amplitude (a), population im- 
balance (b) and the probe intensity (c) versus time when the pump 
amplitude is swept up. The parameters arc the same as those in 
Fig. Q5] The black solid line indicate the result got with method (i) 
while the red dashed line use method (ii). 



Figs. EJa) and (b) are for rj/n = 0.5, the other pa- 
rameters are the same as those in Fig. [31 In this case, 
the probe field intensity is relatively low. As a result, 
the atomic dynamics show little deviation from that de- 
scribed by Eq. ([5]). However when we increase t]/k to 
the value of 1.5, the oscillation period of z (t) experience 
a dramatic enhancement compared to that predicted by 
Eq. ©, as can be seen from Fig. Etc). This critical 
slowing down takes place when the system energy ap- 
proaches a critical value of that defined by an unstable 
saddle point, at which the contour line of energy changes 
its topology from a closed to an open line. This behavior 
is usually associated with bist ability and has also been 
predicted in other systems before [2(1 [27| • 

It is self-evident that when the system is swept across 
the bistable region, it cannot stay in the equilibrium solu- 
tions and non-steady state dynamics will be excited. We 
implement this by ramping up the pump amplitude. The 
population imbalance z then experience a sharp transi- 
tion to the localized state and exhibit periodic oscillation 
with most atoms stay in the deep well, as can be seen in 
Fig. [Sib). If r] is swept in the reverse direction, the value 
of z can take the value of less than 0, which means that 
the system has entered into the extended state, as shown 
in Fig. E^b). Comparing Figs. [B]and Figs. [7[ the abrupt 
transition exhibit a hysteretic behavior, which indicates 
bistability and can be readily observed in experiment. 



i 1 r 




FIG. 7: (Color online) Same as Fig. [6] except that the pump 
amplitude is swept down. 



IV. QUASIPERIODIC LATTICE: NUMERICAL 
RESULTS 

From the previous section we can understand that the 
interplay between nonlinearity and disorder may lead to 
bistability between atomic extended and localized states 
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under appropriate conditions. In the more general case 
that is not an integer, we can only resort to numer- 
ical simulation to investigate the steady-state behavior 
of this system. Here we choose j3 = W5 — l) /2 as the 
inverse of the golden ratio. Then the cavity modes form 
a quasiperiodic optical lattice potential for the atoms. 



responding momentum distribution p (k) broadened, as 
can be seen in Figs. [SJe) and (f). The atomic momen- 
tum distribution in a cavity can be readily observed via 
absorption imaging after the time-of- flight .31]. This can 
serve as a clear evidence that the bistability can really 
taken place between the atomic Bloch extended and AA 
localized states. 
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FIG. 8: Bistability between atomic extended and localized states, 
(a) intracavity photon number for the probe mode |a ss | 2 and (b) 
inverse participation ratio P _1 plotted as function of r//n. Atomic 
density (in log scale) and its momentum distribution for the lower 
bistable branch (c)(d) and upper bistable branch (e)(f) at T)/k = 
1.5. The parameters are the same as those in Fig. [3] 

We focus first on the equilibrium properties of the sys- 
tem by numerically solve Eq. (fT]) with 50 lattice sites. 
The numerical results are shown in Fig. |8l The intra- 
cavity photon number for the probe mode as a func- 
tion of the pumping rate is shown in Fig. [Sfa), which 
exhibits hysteretic behavior typical of dispersive opti- 
cal bistability. We characterize the property of atomic 
wavefunction with the inverse participation ratio = 
/ dx \tjj (x)\ 4 — J2j I c j | 4 ' which is shown in Fig. [SJb). 
This quantity reflects the inverse of the number of the 
lattice sites being occupied by the atoms. Hence a larger 
value of P-i means the atoms are more localized in space. 
The atomic density and momentum distribution of both 
the upper and lower bistable branch at rj/n =1.5 are 
shown in Figs. [Hc)-(f). It is clear that in Figs.EJc) and 
(d) the atomic wavefunction in the lower bistable branch 
display typical properties of bloch extended states with 
almost uniform distribution, while the associated mo- 
mentum distribution p (fc) has well-defined peaks located 
at k = 0, ±2fc^. The atomic density profile in the upper 
bistable branch is exponentially localized, with the cor- 
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FIG. 9: The pump amplitude and the probe intensity versus time 
when the pump amplitude is swept up (left column) and down 
(right column). The parameters are the same as those used in Fig. 



We also investigate the transition dynamics from the 
extended states to localized states and vice versa by 
ramping the pump amplitude up or down, the results 
are shown in Fig.[S] Similar to the dynamics displayed in 
Figs.[B]and[71 when the pumping strength is ramped up to 
exceed a critical value, or ramped down to below a lower 
threshold value, the system will become dynamically un- 
stable and the probe intensity is suddenly increased or 
decreased, typical of dispersive optical bistability. 

The corresponding atomic dynamics are shown in 
Figs. [TU] and QTJ Fig. [TUJ is for the case that the pump- 
ing amplitude is ramped up. We start from an extended 
wavefunction, whose momentum distribution has peaks 
located at k = 0, ±2/cl, as shown in the row marked with 
A. Along with the increase of the pump amplitude, the 
probe intensity gradually increases, the atomic spatial 
distribution becomes more inhomogeneous, and a small 
fraction of atoms will be scattered into momentum states 
p = ±2hk p , ±2h (k p — k^) , • • • due to the beating of the 
two standing-wave cavity mode, as shown in the B row. 
When the pumping amplitude becomes so strong that 
the probe intensity greatly increases and the atomic BEC 
break into fragments, this is an intermediate atomic state 
between extended states and localized states, as shown 
in row C. In this process, the inverse participation ra- 
tio finally jumped to a value around 0.05, which is 
much lower than that would be expected for the present 
system in the localized regime (around 0.6 as indicated 
in Fig. [SJb)). Comparing with the steady-state inverse 
participation ratio as shown in Fig.[S][b), one can see that 
the system can follow the instantaneous steady state (the 
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By contrast, when the pumping amplitude is ramped 
down, the inverse participation ratio gradually de- 
crease to a value approaching zero, as shown in Fig. [TT] 
Contrary to Fig. [TUl it displays the behavior of approxi- 
mately following the instantaneous steady state adiabat- 
ically despite the occurance of bistable transition. This 
is because that the atoms is initially in a localized state, 
whose overlap with the probe mode is small, as can be 
seen from row A. Reducing the probe intensity only lead 
to the weakening of disorder, the atomic wavefunction 
diffuses and the sites neighbouring to the localized sites 
start to be macroscopically occupied, as indicated in row 
B and C, and the atoms becomes more and more delocal- 
ized, roughly following the steady instantaneous steady 
state. 



FIG. 10: Evolution of the inverse participation ratio P when 
the pump amplitude is swept up (the Uppermost figure) . The sub- 
sequent figure show the atomic density (in log scale) (left column) 
and its momentum distribution (right column). The row marked 
with A, B and C correspond to uiRt = 100, 250 and 350, respec- 
tively. 

lower branch) initially, but fails to do so after the pump- 
ing rate reaches the turning point of the lower branch. 
Beyond this point, the probe mode intensity jumps up 
and the atomic wavefunction, instead of forming a lo- 
calized wavepacket as in the case of the steady state, 
becomes fragmented. 
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FIG. 11: Same as Fig. 1101 except that the pump amplitude is 
swept down. Row A, B and C correspond to ujRt = 100, 250 and 
400, respectively. 



V. CONCLUSION 

In summary, we have investigated the system of a non- 
interacting BEC dispersively coupled to two standing- 
wave cavity modes. The equilibrium properties of the 
system are studied under the mean-field approximation. 
In the special case of 1//3 = 2, the cavity modes form a 
lattice of asymmetric double-well which can be well de- 
scribed by a two-mode model. The results obtained from 
the two-mode model and the self-consistent imaginary- 
time propagation of the GP equation show very good 
agreement. Due to the interplay between the nonlinear- 
ity and disorder, bistability between the atomic extended 
and localized states can take place, which is associated 
with the occurance of optical bistability. This is also con- 
firmed via the numerical simulation with an irrational 
in which case the cavity modes form a quasiperiodic 
lattice. 
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